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From a new anti-parallel initial condition using long vortices, three-dimensional turbulence 
forms after two reconnection steps and the formation of at least one set of vortex rings for 
both quantum and incompressible vortices. The long domain allows multiple reconnect ions, 
which enhance vortex stretching rates and the generation of small-scale vortex structures 
within the vortex rings. For the Navier-Stokes vortices, further new features are a profile 
less likely to shed vortex sheets and an improved mapping the direction of the vorticity 
onto the three-dimensional mesh. The vortices evolve via the following steps: First, until 
the first reconnection, dynamics largely consistent with how vortices attract with a possible 
singularity in the Euler equations. Second, vortex reconnection where the symmetry planes 
meet. Third, a series of vortex rings, with the stretching at each following set of reconnections 
leading to the new reconnections and rings. Half of the circulation reconnects into two 
"bridges", leaving two "threads", which the extra stretching transforming the threads into 
spirals wrapped around the bridges. It is argued that these spirals are the source of the 
observed k~ 5 ^ 3 energy spectra and other statistics commonly associated with high Reynolds 
number turbulence. 

PACS numbers: 

To appear in the Procedia IUTAM volume of papers for Understanding Common Aspects of 
Extreme Events in Fluids, Dublin, 2-6 July 2012. as Fully developed hydrodynamic turbulence 

from a chain reaction of reconnection events. 

I. INTRODUCTION 

This paper presents initial comparisons between new Navier-Stokes reconnection calculations 
and recent quantum vortex reconnection calculations. It is also the first of a series of reports on 
hydrodynamic Navier-Stokes and Euler simulations using a new set of anti-parallel vortex initial 
conditions. Figure [I] shows the initial configuration of vorticity used in all of these calculations, 
with vortices that are very long compared with the size of the initial perturbation. For the classical 
calculations isosurfaces of vortices and vortex lines are shown. For the quantum case isosurfaces 
of very low density are shown. 
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Preliminary analysis of these new simulations provides insight into these questions: 

a) Could a simple configuration of vortices generate turbulence promptly from a clearly identi- 
fied series of steps? 

b) Can interacting vortices directly generate a fc -5 / 3 energy spectrum without the usual statis- 
tical cascade arguments? 

c) What is the role of temporal intermittency in determining the higher-order statistics of 
turbulence? 

The conclusion of this paper will be that at least two reconnection steps are needed in the 
interaction between the two vortices before the classical signatures of fully-developed Navier-Stokes 
turbulence appear. One of these signatures is a A: -5 / 3 energy spectrum. The development of -5/3 
spectrum appears to be linked to the generation of stretched spiral vortices [13] , which in turn grow 
out of locally orthogonal vortices created by a series of vortex reconnection events, a configuration 
that has been seen in some isotropic calculations [6l[8] just before they become fully turbulent. 

This paper is organised as follows. First, the new initial condition. This will allow alternative 
numerical methods to reproduce the results given here. Third, the stages in the evolution of the 
vortex structures are described using figures [I] and III, starting from the first reconnection event 



and ending with a turbulent state of vortex rings with swirl. Some statistics for enstrophy growth 
will be provided to give a time-line for the evolution. After the final structure is described, spectra 
and then vorticity moments will be discussed to illustrate the evolution towards a turbulent state. 



II. BACKGROUND: CLASSICAL VERSUS QUANTUM DECAYING TURBULENCE 

Experimentally, a quantum vortex tangle, that is quantum turbulence, is observed to decay 
[El [18] with the same power law as homogeneous, isotropic incompressible turbulence in a periodic 
box and is observed to have a fc -5 / 3 , as does classical turbulence [THE]. Numerically, this can be 
achieved using isotropic initial conditions [19]. However, until recently [10] it has been impossible 
to achieve this type of quantum turbulent state numerically when the simulation starts from only 
two interacting vortices. Which is also true of the classical case. The first goal of the calculations 
reported here was to modify the successful quantum initial condition of very long anti-parallel 
vortices [10J for the Navier-Stokes equations, then to see if this could also evolve into a turbulent 
state. 
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FIG. 1: Comparison of long anti-parallel evolution for Gross-Pitaevskii and Navier-Stokes. The objective 
was to have an initial condition where vortex lines did not separate from the vicinity of the prescribed 
trajectory. In each case a series of reconnect ions results in the formation of a series of vortex rings. The 
Navier-Stokes rings swirl due to the incomplete reconnection of circulation. The steps are discussed in the 
text. 

A. Equations 



The respective classical and quantum equations are as follows. First their respective velocity 
equations ( 1|3 ) are given, which show the similarities between the two dynamical systems. Much 
of the mathematical analysis of these equations is done with the vorticity lj = V x u and the 
wave function ^ = y/pe*^, whose equations are ( 2|4 ) respectively, and do not resemble one another. 
For classical systems either ([!]) or ^ can be integrated in time, for the quantum case only Q is 
integrated. 
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Incompressible Navier- 
Stokes/Euler: 

Gross-Pitaevskii 

^ + (u.V)u=-Vp+ vAu (1) Qp 

dissipation P ( + ( u ' V)u J = -Vp + VE and — + V • (pu) = 

using p _ 1 or V u where u = V0, p = 0.25p 2 , S = a quantum tensor. 



What can be integrated: vorticity what is integrated: ^ = ^pe 
uj = V x u 

0.5V 2 ^ + 0.5^(1 - |^| 



— +(u • V)u? (u? • V)u + z/Au; 
at v v ' n v ^ ^-^-^ 

advection vortex stretching dissipation 

(2) 

The essential difference between the two systems is that the Navier-Stokes equations are dissi- 
pative, while the Gross-Pitaevskii equations are an energy conserving Hamiltonian system. 



B. Initialisation 

For the three-dimensional incompressible equations, while vorticity is calculated by taking 
derivatives of the velocity, one can still initialise with vorticity by inverting its fields to gener- 
ate the initial velocities that are usually integrated in time. However, this is easier said than 
done. 

These issues underlie some of the recent controversy over whether the Euler equations can, or 
cannot, have singularities, as discussed in the issue of Physica D containing the proceedings of the 
2007 Euler 350 meeting, including j2J. The problems with earlier anti-parallel initial conditions, 
going back to |15| . are the appearance of negative spots of circulation in the primary symmetry 
plane, as shown by figure IIB[ left). The consequences of those negative spots discussed in [2] and 



the older conferences proceedings it mentions. Simple fixes such as in figure II Bright) have not 
provided satisfactory corrections. The inspiration for the new initial condition outlined here only 
came about by attempting to apply an inverted classical initial condition to quantum fluids, which 
was then replaced after it developed unphysical instabilities [TO] • The three primary elements of 
the new initial condition are these: 

i) A new profile of the vorticity distribution in the core that is based upon the Rosenhead 
regularisation of a 2D point vortex: 

v(x) = T , ^-^-^ 2 and „Vx,- 



(x - x ) 2 + (y- yo) 2 + a 2 ((x - x ) 2 + (y - yo) 2 + a 2 ) 2 ' 

(5) 



arXiv:1212.3829 



5 



Compare this profile with the Fetter Pade approximate for the density Q. Similar except 
this has been taken to the fourth power. This smooth profile behaves differently than the 
squared-off Gaussian profile of vorticity introduced by [15J and used too many times since 

mmm- 

ii) A new direction algorithm whose goal is to find the nearest position (x 5 , y s , z s ) on the given 
analytic trajectory to a given (x^, yj, Zk) on the three-dimensional grid. The direction of the 
vorticity at the points {x^y^z^) is then given by the tangent of the chosen trajectory at 

iii) Making the vortices very, very long to minimise boundary effects and allow several recon- 
nection events. A new goal is to have a completely local perturbation that asymptotes 
into a straight line. In this way, if the vortex line is lengthened, the perturbation remains 
unchanged. The double cosine path of the past Euler work [2j [9] did not satisfy these 
conditions. In contrast, the following analytic trajectory does. 

x u {x,y,z) = (s x [2/cosh([y/S y } Pt ) - l],l,o) . (6) 

The following lengths are used: 8 X — —1.6 and 8 y — 1.25. pt = 1.8 was chosen in order to 
localise the perturbation near the y = symmetry plane. 

• The different physics in the quantum case dictates how to choose the density and phase 
around the p = vortex cores. The density profile p(r) must go from p(0) = to p — >> 1 as 
r — » oo. One example is this low-order Pade solution [3]: 



^(r) =r/Vr 2 + 2 or p(r) = r 2 /(r 2 + 2) . (7) 
Points i) and ii) resolve the problem with the negative spots, leading to the final profile in figure 



II B This has been used for anti-parallel vertical vortices in a stratified fluid, the anti-parallel 
unstratifled vortices mentioned here, and anti-parallel Euler vortices. Unphysical initial instabilities 
due to small-scale inbalances have been suppressed, without the extra massaging or squeezing used 
before 019], resulting in calculations where a cleaner and stronger larger-scale instability forms. Is 
this sufficient to lead to a transition to sustained turbulence from vortex interactions where none 
had been seen in earlier work? 

The two simulations are used here say it is. Removing the symmetries, the full domains are 
L x x L y x L z = Atv x 327T x Att and L x x L y x L z = 4tv x 32tv x 8tv, with x, y, z values going 
as approximately [0 : 12, —50 : 50, —6 : 6] and x, y, z going as [0 : 12, —50 : 50, —12 : 12]. The 
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Melander-like jell=0 jtype=0 jprof=0 T=0 




FIG. 2: Left: Noise generated by one of the original anti-parallel initial conditions [15]. Right: From 
[2]. As indicated by the position of the min(co> y ). As in [9], r±_ is used along with a squared-off Gaussian 
profile, which is then lengthened in x, in this case by using an anisotropic filter in k x . In [9] the profile was 
lengthened by adding a localised du x /dz shear. 

L z — 87T, v — 0.001 calculation is used to show the full spatial development of the structures and 
the L z — 47T, v — 0.0005 case is used to show that a full Kolmogorov A: -5 / 3 inertial subrange 
with very strong — S u fluctuations develops. The important timescales are the same in the two 
calculations. 

The insets show the upper/left quarter domain near the y = perturbation plane. The t = 
inset shows that the circular cross-section, taken perpendicular to the direction of vorticity Co in its 
centre, has constant width along its entire length and the t — 16 inset shows that large fluctuations 
of vorticity about zero seen previously [2j [7] are now avoided. 

All of the initial conditions to be discussed then apply a Fourier filter of the form 

with 84 = 0.01. This initialisation was done on modest meshes, typically 128 3 , before those fields 
are remeshed onto the meshes of true interest. 

III. INSTABILITY STEPS 

Initial instability. The primary inviscid linear instability on anti-parallel vortices was orig- 
inally identified by [3J. This quickly becomes nonlinear with the perturbation in anti-parallel 
vortices in figure [I] evolving into the structures GP t = 1.5 and NS t = 16 structures that are 
beginning to reconnect at the point of closest approach while generating waves away from those 
points. The insets at t = and NS t = 16 show the perspective that has appeared in a number of 
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FIG. 3: Left: x — y cross-section of uo at t = for y > through the position of maximum vorticity, which 
is on y = near x = 7.25. by the continuous crosses of uj y = at these boundaries. Right: Circulation 
in the two primary half-symmetry planes. Their sum remains constant through the reconnection process 
which results in nearly half of the original circulation being reconnected. 

earlier papers, in particular [9J, with a cut through the symmetry plane. This cut shows familiar 
features including the leading heads, remnants of the original vortices, and the following tails, 
flattened vortex sheets. These sheets are why the extra resolution is needed in z. 



After t = 16, as the reconnection continues, the twists turns into kinks, shown in figure III 
These twists and kinks do not form if periodic perturbations are used [15J, which the long domain, 
plus the localised perturbations, promote. 

Ongoing analysis of the curvature for t < 16 has found a complicated structure of peaks and 
dips along the vortices with the curvature on the y = perturbation plane remaining small for all 
times. These peaks and dips are within the isosurface shown in t = 16 frame and are not directly 
connected to the bulges to larger \z\ near \y\ = 5. These bulges are more directly related to the 
requirement that any volume that has been pulled away from the y — perturbation plane by 
vortex stretching must go somewhere, and in this case it goes into the bulges. Another feature of 
the bulges is how they bend back towards the z = dividing plane, overshooting the original z 
positions of the vortices. 

By t = 32, the first reconnection event has completed, with the heads reconnecting across the 
z = plane to form bridges and kinks that develop from the inviscid twist generated before t = 16. 

Roll- up and twist As the first reconnection is finished at t = 32, the tails roll- up to form 
new tubes on the y = perturbation plane. These new tubes are commonly known as braids or 
threads, with one indicted by the red vortex line coming out of the y = perturbation plane in the 
t = 96 frame. To be able to observe the complicated interactions at these two times, two vorticity 
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FIG. 4: t — 32: Evolved from t — 16 when reconnection began. Left: Two isosurfaces. Outer uj = .63 
transparent surface in turquoise for comparison with t = 16 ang later times, and an inner uj = .89 surface 
to capture the finer vorticity details. Right: The u = .55 isosurface is shown to emphasize how the tails 
have started to roll-up, seen on the right of the cut through the y = perturbation plane. How the roll-up 
connects to the original vortices is in red while the reconnected parts of the vortices are blue. They are 
following each other into the main structure, but are not yet twisting around one another. 



isosurfaces are used: \ou\ = 0.89 and \ou\ = 0.63, shown only for the y > half-domain so that the 
y = 0, x — y cut is visible. 

The two important features highlighted by the |cj| = 0.89 isosurface at t = 32 are the roll- up of 
the tails and the initial twisting between the bridge and tails. At t = 32, the new threads tuck behind 
the reconnected bridge at y « 2 and z = ±1 with the two sets of vortices are aligned orthogonally, 
similar to what has been seen after vortex sheets first reconnect in simulations initialised with 
only a few small wavenumber Fourier modes [U [6]. Beyond the orthogonal crossing point, the 
unreconnected threads -htails begin to wrap around the reconnected bridges, near y w 3, generating 
a twisted structure that has not been unraveled further. 

The lower \u\ = 0.63 surface shows the relationship between these intense regions and the new 
path for the envelope of the original vortices. This includes the bulge that goes to \z\ > 4 near y = 5. 
Note the new intensification of vorticity inside the \u\ = 0.63 isosurface for 5 < \y\ < 10. This 
growth is being driven by a new anti-parallel interaction that is developing around 8 < \y\ < 12. 
This new source of vortex stretching plays a role in the next stage illustrated by the frame at 
t = 96. 

Timescales In order to interpret the time evolution shown by the series of frames in figure [TJ 
diagnostics are needed for identifying timescales. One comes from normalising the production of 
enstrophy Z into the following skewness: 

dZ/dt 



/9 where [S]ij — - ( + ) an d c i s an isotropy coefficient . (9) 
/j i z y CJXj cjx<i 
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FIG. 5: Vorticity and energy time evolution for the v = 0.0005, L z = 2ir calculation. Left: Decay of kinetic 
energy (lower /black), growth until t = 320 (then decay) of enstrophy and velocity derivative skewness: 
—S u . Several intermediate timescales set using the growth then decay of the enstrophy Z and —S u , its 
normalised production. First timescale: t\ ~ 48, showing the largest — S u and representing the completion 
of the first reconnection, which began about t c ~ 14. t<i ~ 110 is completion of the second reconnection and 
detachment of the first vortex rings. The —S u « 2 overshoots the steady-state value of about — S u = 0.5 to 
0.7 by an unprecedented margin compared to earlier work. But is not due to singular growth, which Euler 
calculations say should end after t ~ 14. Right: Temporal variation in the higher-order vorticity moments, 
normalised and scaled in the manner suggested by [5] . By using this scaling on this modestly large Reynolds 
number Navier- Stokes calculation, one can see that each higher-order moment is bounded from above by 
the lower-order moments. 



It is known from experiments and simulations that at large Reynolds numbers — S u — >> 0.5 — 0.7, 
so the growth and relaxation of — S u to this asymptotic value can used to set the timescales of the 
calculation. 

Figure [5] shows the most significant timescales appearing in the global diagnostics. On the left 
is the evolution of the global enstrophy Z, its normalized production — S u Q and the decay of 



kinetic energy. On the right is the evolution of the rescaled vorticity moments D m (10). The 



D m are characterised by a hierarchy where the lower-order D m bound the higher-order D m for all 
times and order. This hierarchy, including comparisons of these curves to their counterparts from 
high-Reynolds number homogeneous, isotropic simulations, is being investigated further. What is 
important for the evolution described here are the following. First, the growth of the Z, — S u and 
the D m are nearly identical up until t w 16. This is when the first reconnection begins and is a 
candidate time for when there could be a singularity in the v = Euler limit. 

The next important time is t « 48 when —S u has its first and largest peak for both cases and 
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there is a minor peak in Z. Isosurfaces at t — 48 are almost the same as those at t = 32, with a 
little additional twisting up the reconnected bridges. —S u decreases rapidly for t > 48, reaching 
a minimum at t w 90, which is when the second reconnection is beginning and is discussed next. 
The additional peaks in — S u can all be associated with the additional reconnection events being 
discussed. These features are largely mirrored by the D m in the right frame of figure [5| including 

Spirals By t = 32, the reconnected bridges have begun to separate in y. This is because the 
curvature underlying the Biot-Savart interaction reverses its direction in the reconnected vortices, 
so the interaction goes from attraction in z for t < 16 before reconnection, to repulsion in y for 
t > 32 after reconnection. As the t > 32 vortices separate, the reconnected bridges wind the 
unreconnected threads about themselves, illustrated at t = 96 by how the red thread coming from 
the y = plane snakes around the blue line coming from the z = bridge. For large y, the red and 
blue lines meet and joint into the original vortex. 

The other major feature to note in the t = 96 frame is that the new anti-parallel attraction 
seen for 8 < y < 12 at earlier times has now moved out to y w 17 and has progressed to the stage 
where the vortices are starting to touch and reconnect again. The inner isosurface shows flattening, 
curvature, and preliminary roll-up similar to what forms at y = at t = 16, just before the first 
reconnection. Furthermore, for the period just before and after t = 96, the normalised enstrophy 
production and the higher-order D m in figure [5] all grow significantly, as they do around t — 16, 
the time of the first reconnection. 

After the reconnection that starts at t = 96 has been completed, we are left with a disconnected 
vortex ring with strong spiraling patterns snaking completely around it, as shown at t = 176. 
By t = 176, a third reconnection site is developing near y « 35, which leads to a second vortex 
ring disconnecting, with new spirals forming on its low y side. A second ring with spirals then 
disconnects, as shown at t — 256. Given a large enough domain in y, this process would probably 
continue to produce new spiral rings, as in the quantum vortex case [10J. It also has similarities 
to experimental observations such as [12] where thin threads connecting a chain of vortex rings are 
seen. 

IV. SPECTRA AND MOMENTS 

The initial high-wavenumber spectra are given by the hyperviscous filter, which quickly develop 
into power laws that are steeper than A: -4 , as seen in Euler calculations [9J. The subsequent spectral 
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FIG. 6: Kinetic energy spectra: Domain: 27r(2,8, 1) v = 0.0005. Uncompensated spectra shows the pro- 
gression from a k~ 3 spectrum at short times to the formation of a k~ 5 ^ 3 at later times. As the large-scale, 
low wavenumber part of the spectrum begins to decay, Kolmogorov similarity (using the time dependent 
dissipation e(t) = vZ) sees convergence starting at t — 184. The estimated Kolmogorov coefficient Ko is 
between Kq=1.05 and 1.1. This is lower than what is obtained from large forced numerical simulations, 
which typically give Ko « 1.5, perhaps because not all of our global enstrophy Z (and therefore dissipation 
e = vZ) is participating in the turbulent dynamics. 

evolution is shown in figure |6j starting at t = 44. Eventually, starting at t = 184, a A: -5 / 3 begins 
to form, lengthen, and collapse according to Kolmogorov similarity. k x and k z spectra have also 
been determined. They do not have the same extent as the k y direction because the physical space 
domains in those directions are smaller, but their inertial subranges are consistent in magnitude 
and extent with those in the k y direction. 

While spectra can provide some measure of whether a fluid is turbulent or not, they provide little 
information about the statistical nature of physical space velocity fluctuations and intermittency. 
If the experimental conditions are relatively steady for long periods of time, there are a number of 
alternative physical space statistics that can be used to characterise the intermittency of turbulent 
flows. However, determining many of these numerically is restricted by how rapidly they converge 
to steady statistics. In the analysis of high Reynolds number, very large numerical simulation, what 
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[20] have demonstrated is their ratios converge, even if their values don't. As recently pointed out 
[11], this can be explained by the mathematics in [5j. 

What insight can these reconnection calculations provide? They can if the rise and fall of 
these intermittent statistics can be tied to the reconnection events and the growth of enstrophy 
production. 

First, is intermittency growing sufficiently rapidly as the Reynolds number increases for these 
reconnection simulations to be useful? The strong growth in the skewness fluctuations in figure 
|5^1eft) suggests that this is the case as the two initial peaks near t = 40 and t = 105 double when 
the Reynolds number is doubled. Even higher growth is expected when this is redone with an 
L z — Atv domain. 

Which brings us to the rescaling of the higher-order vorticity moments as a function of time 
in the manner recently suggested by [5], given in figure fright). This fi gure shows a hierarchy of 
higher-order vorticity moments, higher than m = 1 (the enstrophy), which would off the upper 
scale at the top. The normalisation is based upon the frequency scale w — v jL 2 z and a power law 
designed to compare the degree of intermittency for each moment and goes as 

D m = (wH m )^ &m where a m = 4m/(4m - 3) , (10) 

where the H m = (J cj^dV) 1 / 2171 are the standard vorticity moments. With this scaling, each of 
the higher-order moments is bounded from above by the lower-order moments. Furthermore, as 
the order increases, the moments become bunched in regions of strong growth. The ordering is 
fostered by the choice of frequency scale and the bunching is fostered by the choice of power law. 
The maximum of vorticity, the m = oc norm, follows the same scaling and appears just below the 
m = 5 curve. 

There is no rigorous proof that this ordering should be obeyed, but this observation should 
be significant for the question of whether the Navier-Stokes equations have singularities or not. 
First, could it be shown that for any sufficiently smooth Navier-Stokes initial condition that these 
moments relax quickly into this ordering, and that this ordering is then maintained? If so, then 
we would know that all the higher-order vorticity norms will eventually be bounded from above by 
the lower-oder vorticity norms. Which brings us to the lowest order, m = 1 norm, the enstrophy. 
Implying that if the growth of enstrophy could be bounded for all times, then all higher-order 
moments of vorticity could be bounded, including Hqq. 
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V. SUMMARY 

The case is made in this paper that if the initial vortices are well-balanced and the domain is 
large enough for two independent reconnection steps to form, then the flow wil rapidly break down 
into three-dimensional, fully- developed turbulence. 

The initial configuration used to demonstrate this consisted of very long, anti-parallel vortices. 
It is argued that the long domain is needed in order to accommodate the multiple reconnections, 
which enhance vortex stretching rates and the generation of small-scale vortex structures within 
the vortex rings. In addition to making the initial vortices very long, several new features are 
applied during initialisation that are designed to make it less likely that the flow will generate 
vortex sheets. These include a new balanced profile and an improved way of mapping the direction 
of the vorticity onto the three-dimensional mesh. 

To get to the turbulent state, the vortices progress through the following steps: First, until 
the first reconnection, the vortex dynamics are largely consistent with existing work on possible 
singularities of the Euler equations. This includes a head and tail structure with the maximum 
of vorticity between these them. A feature not clearly noted before is how the vortices twist back 
upon themselves towards secondary points along the original anti-parallel vortices. This nonlocal 
amplification of the anti-parallel interaction provides extra stretching along the entire vortices, 
including the primary symmetry plane. 

Jumping to the last figure of higher-order moments, note the convergent growth of these mo- 
ments leading up to t = 16. This is related to singular Euler dynamics that is being discussed 
elsewhere. 

The next step is reconnection on the primary symmetry plane. Only the half associated with the 
head reconnects into two "bridges" , leaving the tails to roll- up into two "threads" . Together these 
form a locally orthogonal configuration with steep energy spectra and minimal energy dissipation. 

The third step starts with orthogonal configurations and is modified, compared to past simu- 
lations, by the stretching induced by the anti-parallel interactions at the secondary points. Due 
to this secondary stretching, as the threads begin to wrap around the newly reconnected bridges, 
they are also pulled along the bridges to form spirals. 

Eventually there are two more reconnections at secondary points, which result in vortex rings 
separating from the original vortex lines. It is only at this stage, while the spectra are becoming 
less steep, that significant dissipation develops, first at higher wavenumbers. 

The fourth step will be defined as an extended period with multiple reconnections, forming 
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multiple vortex rings, during which a nascent -5/3 extends to low wavenumbers and additional 
statistics achieve values consistent with fully- developed turbulence. Each reconnection is at the 
secondary interaction positions induced by the previous anti-parallel interaction and reconnection. 
And each pair of reconnections forms a vortex ring. 

In the calculations reported here, there are three distinct reconnection events from which two 
sets of independent vortex rings form on either side of the symmetry plane. For the quantum 
vortex calculations in figure [TJ it is seen that there is no end to the number of rings that can form, 
so long as the domain is large enough. For each set of reconnections, additional stretched spirals 
form as the A: 5 / 3 energy spectrum appears. The final vortex structure is a chain of vortex rings 
connected by thin threads, as seen in some experiments. 

The final state satisfies all the usual properties of fully- developed turbulence, usually associated 
with homogeneous and isotropic flows, neither of which is true for these simulations. There is 
a Kolmogorov A: 5 / 3 subrange with roughly the correct constant. And after highly intermittent 
variations, the velocity derivative skewness — S u , or normalised enstrophy production settles to 
values of about 0.5 to 0.6, which are the values characteristic of fully-developed turbulence. It 
is then shown that the intermittency of — S u can be extended to intermittency in higher-order 
moments of vorticity. Moments that we hope will be determined numerically in some upcoming 
cutting-edge simulations of forced, homogeneous, isotropic turbulence similar to that used by [20] . 
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